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Abstract 

Applying asymptotic expansions at threshold, we compute the two-loop QCD cor- 
rection to the short-distance coefficient that governs the leptonic decay "0 l^l^ 
of a S-wave quarkonium state and discuss its impact on the relation between the 
quarkonium non-relativistic wave function at the origin and the quarkonium decay 
constant in full QCD. 
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Quarkonium decays played an important role in establishing Quantum Chromodynamics 
(QCD) as a weakly interacting theory at short distances. Calculations of heavy quarko- 
nium decays usually proceed under the assumption that the heavy quark-antiquark 
bound state is non-relativistic and that the decay amplitude factorizes into the bound 
state wave function at the origin and a short- distance quark-antiquark annihilation am- 
plitude. One can then explain the small width of the S-wave spin-triplet charmonium 
state J/ip, because it can decay only through electro-magnetic annihilation or annihila- 
tion into at least three gluons Today's understanding of quarkonium bound states 
has refined this picture and allows us to calculate relativistic corrections systematically 
at the expense of introducing further non-perturbative parameters that characterize the 
bound state. Such calculations can be done most transparently in the framework of 
a non-relativistic effective field theory (NRQCD) [Q,^ that implements the factoriza- 
tion of contributions to the (partial) decay widths from different length scales. Besides 
potentially large relativistic corrections, the size of radiative corrections to the quark- 
antiquark annihilation amplitudes has always been a matter of concern. The one-loop 
radiative corrections to the decays ip —>■ l^l^ (where / = e,/i) 0], ^/^ —flight hadrons, 
if) — > 7+ light hadrons 0] are large and question the practicability of factorization for 
charmonium and, perhaps, even bottomonium. [Here and in the following we use if) as 
a label for any S'-wave spin-triplet state, i.e. J/ip or ip' for charmonium and T{nS) for 
bottomonium.] 

In this Letter we address this question and report on the calculation of two-loop 
short- distance corrections to the leptonic decay ip l^l^ . This is the first (and also 
'easiest') two-loop matching calculation for quarkonium decays and it can also be used to 
connect the decay constant of the meson defined in QCD with the non-relativistic wave 
function at the origin, which appears in NRQCD and potential models. In turn, this 
wave function is an important input parameter for the prediction of other quarkonium 
decays and also quarkonium production cross sections. 

We recall that if) decays leptonically through interaction with the electro-magnetic 
current. The partial decay rate, neglecting the tiny lepton masses, is exactly given by 

where is the mass of ip, the fine structure constant and eg the electric charge 
of the heavy quark in units of the electron charge. The ip decay constant is defined 
through the following matrix element of the electro-magnetic current: 



[e^{p) is the ip polarization vector and p the ip momentum.] The leptonic decay rates are 
known experimentally [0. For J/ip we find fj/^ = (405 ± 15) MeV. 

The decay constant parametrizes the strong interaction effects and contains long- 
and short- distance contributions. For quarkonium the short-distance scale is l/M^ and 
the long-distance (bound state) scales are l/{M^v) and l/(M^t>^), where v is the (small) 
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characteristic velocity of the ip^s quark constituents. The short- distance contributions 
can be isolated, and calculated in perturbation theory, by matching the vector current in 
QCD onto a series of operators in NRQCD. Up to corrections of order v^, the matching 
relation is given by 



Co «„^)(^|^W|0)(/i) 



Qm3 



(3) 



where (not to be confused with the ip meson) and x denote non-relativistic two-spinors, 
D the spatial covariant derivative and tuq the heavy quark mass. [More details on 
notation and NRQCD can be found in Ref. 0. Note, however, that we use a relativistic 
normalization of states also for the matrix elements in NRQCD.] The matrix elements on 
the right-hand side are defined in the ip rest frame and A{p) is the matrix that performs 
the Lorentz boost into this frame. The matching coefficients Cq and Ci are expressed 
as series in the strong coupling as and account for the short-distance QCD corrections. 
The matching coefficients and matrix elements in NRQCD individually depend on the 
factorization scale /i. Ci is defined (as in Ref. [0) such that Ci = 1 + 0{as) and 0] 

C„ («,. ^) = 1 - ?£l^ + {^)\-, (4) 

where Cp = {N^ — l)/{2Nc) and Nc = 3 is the number of colours. We now discuss the 
calculation of the two- loop matching coefficient C2. 

Since the matching coefficient contains only short- distance effects, it can be obtained 
by replacing the quarkonium state ip on both sides of Eq. (^ by a free quark-antiquark 
pair of on-shell quarks at small relative velocity. In terms of this on-shell matrix element, 
the matching equation can be rewritten as 

^2, QCD Tqcd = ,NRQCD Zj rNRQCD + 0{v'), (5) 

where Z2 are the on-shell wave function renormalization constants in QCD and NRQCD 
and r the amputated, bare electro-magnetic annihilation vertices in QCD and NRQCD. 
The Feynman diagrams for Fqcd at two loops are shown in Fig. [l|. Since the current 
J = tp^aix need not be conserved in NRQCD, we allowed for its renormalization, Jbare = 
ZjJrcn, on the right-hand side. We then obtain Cq by calculating all other quantities in 
Eq. (^ in dimensional regularization and using the modified minimal subtraction scheme 
(MS scheme) [|. 

The matching calculation is considerably simplified, if one uses the threshold ex- 
pansion of Ref. to compute Fqcd directly as an expansion in v^. The threshold 
expansion is obtained by writing down contributions corresponding to hard (/ ~ fnq), 
soft (/ ~ rnqv), potential (/q ~ rnqv^, U ~ mqv) and ultrasoft {I ~ rnqv^) regions. The 
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Figure 1: Diagrams that contribute to Tqcd- Symmetric diagrams exist for -02.3,5- The last 
diagram summarizes vacuum polarization contributions from massless fermions (Dq), gluons 
(-D7), ghosts {Dg) and the massive fermion with mass niQ {Dg). 

contributions from soft, potential and ultrasoft loop momenta can all be identified with 
diagrams in NRQCD that appear in the calculation of rNRQCD- Hence, they drop out of 
the matching relation Eq. @ and it suffices to compute the contribution to the threshold 
expansion of the diagrams in Fig. |I], where all loop momenta are hard. [The threshold 
expansion is not only convenient; it also provides an implicit definition of NRQCD in di- 
mensional regularization. This is necessary, because dimensionally regularized NRQCD 
is not given by the dimensionally regularized Feynman integrals constructed from the 
vertices and propagators of NRQCD. In order to avoid that the cut-off for the effective 
theory be treated as larger than rriQ, dimensionally regularized NRQCD has to be sup- 
plemented by a prescription for expanding the Feynman integrands. This prescription is 
provided to all orders as part of the diagrammatic threshold expansion method 0.] 

We briefly describe the calculation of the hard contributions to Fqcd- [Details of the 
calculational method and a solution to the recurrence algorithm for the two-loop integrals 
will be given in a long write-up of this Letter.] The spinor structure of the on-shell matrix 
element in QCD is conventionally parametrized by two form factors, Fi and F2, of which 
only the combination Fi + F^ is required here. Since terms of order f ^ are not needed 
to determine Cq [see Eq. (^], we may set the relative momentum to zero and compute 
the form factors directly at threshold. The form factors have Coulomb singularities at 
threshold and diverge as However, these singularities appear only in the soft, 

potential and ultrasoft contributions, and the hard contribution is well-deflned directly 
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Table 1: Coefficient of (as/Tr)^ (e'E / (^4^^ ^2^y2e f^j, ^j^g ^lard contribut ion to the diagrams 
of Fig. H] evaluated at threshold = 4m^. 02^3^5 include a factor of 2 to account for the 
corresponding symmetric diagrams. For SU{3) the colour factors are Cf = 4/3, Ca = 3, 
Tp = 1/2. The number of light (massless) quark flavours is denoted by n/. 



at threshold in dimensional regularization. The loop integrals simplify considerably, once 
the relative momentum is set to zero, since they then depend only on a single scale. We 
then project on the form factor Fi + F2 and reduce all integrals to integrals without 
numerators. These integrals can be further reduced to 'simple' integrals and two non- 
trivial two-loop integrals by means of recurrence relations derived from integration by 
parts in the loop momenta [|l^. The solution to the recurrence relations is the difficult 
part of the calculation. The remaining non-trivial two-loop integrals can be calculated 
explicitly using standard Feynman parameters. The results obtained for the diagrams 
of Fig. |1| are listed in Table ^ 

After summing all the diagrams, multiplying by the two-loop QCD on-shell wave 
function renormalization constant |rT|, and performing (one-loop) coupling and mass 



renormalization the result still contains poles in e = (4 — (i)/2 (where d is the space- 
time dimension). Since the wave function renormalization constant 2^2,nrqcd in NRQCD 
equals 1 up to higher-order terms in v^, not needed here, these poles are attributed to 
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an anomalous dimension of J, which first arises at the two-loop order. As a consequence 
the matrix element ('?/'|'?/'Vjx|0)(/i) is factorization scale-dependent. We define it in the 
MS scheme and obtain the anomalous dimension for the NRQCD vector current J: 



7j 



dlnZj 
din /I 



+ 0{a 



(6) 



The scale dependence is compensated by the scale dependence of the two-loop matching 
coefficient C2(mQ//i) in Eq. (H). Separating the different colour group factors, the final 
result for C2{mQ/^) in the MS scheme is: 



Ciimq/fj.) = Cp C2,A + CfCa C2,na + CFTFnfC2,L + CfTf C2,h, 
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Here C(3) = 1.202 . . . and we have taken all fermions with masses less than ttlq as 
massless, which is a good approximation even for mg = rrih, the bottom quark mass, in 
which case we neglect rric, the charm quark mass. The coefficient C2,l proportional to Uf, 
the number of light fermions, has been recently obtained by Braaten and Chen ^ . Note 
also that the Cp-teim of the form factors Fi 2 close to threshold has been calculated by 
Hoang using the absorptive parts of the form factors obtained in Ref. [0. Because 
the result in Ref. |]12| contains hard and soft (potential, ultrasoft) contributions, it is not 
possible to extract the matching coefficient C2,a from Ref. [|T2|. [The structures in C2,a 
that cannot arise from the small loop momentum regions agree with Hoang's result.] 

The size of the two-loop correction to Co [Eq. (^] given by Eqs. (0)-(|ll|) is enormous. 
We define the (scale-dependent) non-relativistic decay constant as 

iijl^P^a.xmfi) = (-0/f^^(/x)M^e*, (12) 

in analogy with Eq. (^. The non-relativistic decay constant is related to the ip wave 
function at the origin by {f^^Y = 12 |^(0)p. Using Eq. (g), we obtain 



1 - 



jas(mQ) 
37r 



- (44.55 - OAlUf) [ 



+ 



(13) 



With as{mc) ~ 0.35 and as{mb) ~ 0.21 the second-order correction exceeds the first- 
order correction even for the bottomonium states. For charmonium, the second-order 
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term is almost twice as large as the already sizeable first-order correction. [Note that 
the BLM estimate of the two-loop correction is far off the exact two-loop result.] 
Perturbative matching at the scale fi = mq does not seem to work. Can the factorization 
of short- and long-distance effects still be useful? 

A novel, and perhaps unexpected, aspect at the two-loop level is the factorization 
scale dependence of the non-relativistic decay constant and, hence, the quarkonium wave 
function at the origin. The scale dependence is large, especially due to the non-abelian 
term in Eq. (|^). This scale dependence of the wave function indicates the limitation of 
the non-relativistic potential model approach already at leading order in f ^, since the 
wave functions obtained from solving the Schrodinger equation are scale- independent. 
Despite this shortcoming it may be argued that the wave function/non-relativistic decay 
constant obtained from potential models corresponds - if to anything - to the wave 
function evaluated at a scale typical for the bound state and not mq. This point of 
view is also evident if the wave function is computed non-perturbatively using lattice 
NRQCD, in which case the ultraviolet cut-off/factorization scale is also much smaller 
than mq. We consider /i = 1 GeV as an adequate bound state scale for bottomonium 
and charmonium, as the applicability of perturbation theory prevents us from taking yet 
smaller scales. We then find, for bottomonium {rrih = 5 GeV, Uf = 4): 

fr,.s> = (l - ^ - 1.74 {^^^ \...) fZsfi GeV) « 0.81 f^^fi GeV). 

(14) 

The numerical factor 0.81 is stable against variations of the scale of the coupling at 
fixed factorization scale 1 GeV. At this factorization scale the second-order correction is 
numerically insignificant. Although we do not know whether the three-loop correction 
would also be small at the low factorization scale, we tend to consider Eq. (|l^) as an 
accurate prediction. This prediction could be tested, if /T(n5)(l GeV) were accurately 
known, for example from NRQCD lattice simulations. The large scale dependence raises 
the theoretical question (the answer to which we postpone to the long write-up) as to 
whether it is necessary to resum the logarithms ln(mQ//i^) to all orders. 

The situation is less favourable for charmonium states. Since rric ~ 1.5 GeV, the 
size of the second-order correction is altered little for scales ^ > 1 GeV. We have not 
succeeded to find a trustworthy interpretation of Eq. (|^) and conclude that the fac- 
torization of short-distance and long-distance effects may not be useful in practice for 
charmonium. Since the leptonic decay is the simplest conceivable decay, this puts into 
question the possibility to obtain universal relations between various charmonium decays 
and production processes through the use of NRQCD [^Q- This pessimistic conclu- 
sion may be biased by our use of the MS factorization scheme. It is conceivable that 
other factorization schemes or relations between physical observables, from which the 
wave function at the origin is eliminated, exhibit better convergence properties of their 
perturbative series. A definitive conclusion on this issue can be obtained only once a 
second quarkonium decay is computed to two-loop order. 
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While this paper was being written, Czarnecki and Melnikov also considered the 
two- loop matching of the electro-magnetic current, also using the asymptotic expansion 
method of 0. After correction of a trivial normalization error (A. Czarnecki and K. 
Melnikov, private communication), their result agrees with ours. 
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